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Decoherence mechanisms in crystals of weakly anisotropic 
magnetic molecules, such as V15, are studied. We show that 
an important decohering factor is the rapid thermal fluctua- 
tion of dipolar interactions between magnetic molecules. A 
model is proposed to describe the influence of this source 
of decoherence. Based on the exact solution of this model, 
we show that at relatively high temperatures, about 0.5 K, 
the quantum coherence in a V15 molecule is not suppressed, 
and, in principle, can be detected experimentally. Therefore, 
these molecules may be suitable prototype systems for study 
of physical processes taking place in quantum computers. 
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A new class of magnetic compounds, molecular mag- 
nets [Q, has been attracting much attention. Each 
molecule of such a compound is a nanomagnetic en- 
tity with a large spin (or, in the antiferromagnet case, 
large staggered magnetization). The interaction between 
different molecules, being of the dipole-dipole type, is 
very small, so that the corresponding crystal is an ar- 
rangement of identical weakly interacting nanomagnets. 
Molecular magnets are ideal objects to study phenomena 
of great scientific importance for mesoscopic physics, such 
as spin relaxation in nanomagnets, quantum tunneling of 
magnetization, topological quantum phase interference, 
quantum coherence, etc J2|. Low-spin weakly anisotropic 
compounds, like V15 [g|4], demonstrate well-pronounced 
quantum properties, such as significant tunneling split- 
ting of low- lying spin states. As we show here, they are 
attractive prototype systems to study mesoscopic quan- 
tum coherence and physical processes which destroy it. 
Besides fundamental science, these studies are important 
also for the implementation of quantum computation M . 

At present, for strongly anisotropic high-spin magnetic 
molecules such as Mni2 and Fes, different kinds of de- 
cohering interactions have been studied p,[7| and their 
interplay with quantum properties at low temperatures 
(below f.5-2 K) is well understood. A general conclu- 
sion M about strongly anisotropic systems is that the 
dissipative environment, especially the bath of nuclear 
spins, rapidly destroys coherence even at very low tem- 
peratures, and only incoherent tunneling survives. 

Decoherence in weakly anisotropic magnetic molecules 
has not yet received much study, and such a study is 
the main purpose of the present paper. We analyze var- 
ious sources of decoherence for such molecular magnets 
as V15, and show that in the temperature range 0.2 K 
- 0.5 K the decoherence is governed by rapidly fluctuat- 



ing dipole-dipole fields produced by thermally activated 
molecules. This mechanism in molecular magnets has not 
been considered before; estimates show that in strongly 
anisotropic magnets like Mni2 or Feg this effect is small. 
Based on an exactly solvable model, we demonstrate that 
even at temperatures as high as 0.5 K, the quantum co- 
herence in V15 molecules is remarkably robust, and, in 
principle, can be detected experimentally. Therefore, the 
V15 molecular magnet is a promising candidate for the 
study of quantum coherence and may be a useful pro- 
totype system for the investigation of physical processes 
taking place in quantum computers. 

The magnetic subsystem of the molecule 
K 6 [Vi5As 6 042(H20)]-8H 2 0, (denoted for brevity as V15) 
consists of fifteen V 4+ ions with the spin 1/2 (see Fig. nl). 
The ions form two nonplanar hexagons (with total spin 
equal to zero) and a triangle sandwiched between them. 
Exchange interactions between ions are reasonably large 
(from 30 K to 800 K), but, due to the strong spin frus- 
tration present in the molecule, the couplings of the cen- 
tral triangle spins with the hexagons cancel each other 
(see Fig. Q). The hexagon spins form a rather stiff anti- 
ferromagnetic structure, and the low-energy part of the 
spectrum is defined by only three weakly coupled spin 
1/2 ions belonging to the central triangle. An effective 
exchange coupling between the triangle spins J ~ 2.5 K 
is present. Thus, the ground-state term consisting of two 
doublets with S — 1/2 is separated from the low-lying 
excited term S = 3/2 by the distance Ai = 3J/2 ~ 3.8 
K. Experimental results [Q suggest that within the two 
ground-state doublets, the states \S = 1/2, S z = +1/2) 
and \S — 1/2, S z = —1/2) are mixed (a small anisotropic 
interaction may be responsible, but it is not of concern 
for the arguments presented), so that tunneling between 
these levels occurs and the fourfold degeneracy of the 
ground state is partly lifted (it cannot be lifted com- 
pletely because of Kramers' theorem: in the absence of 
an external field all levels are doubly degenerate). The 
coherent tunneling leads to a splitting Ao — 0.2 K pj 
between the two pairs of Kramers-degenerate levels. The 
aim of this paper is to study the decoherent influence 
of the environment upon this tunneling, i.e. the deco- 
herence between the states |5 = 1/2, S z = +1/2) and 
|£= 1/2,5, = -1/2). 

First, we consider decoherence caused by the spin- 
lattice relaxation. The rate of the relaxation due to direct 
one-phonon processes can be estimated MM as 
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where Ao is the tunneling splitting of the ground state 
doublets, T is the temperature, v « 2800 m/s Q is the 
sound velocity in the molecular crystal, M « 2.3 • 10 3 
a.m.u. is the mass of the molecule, = (67T 2 i; 3 /r2o) 1 ^ 3 — 
70 K is the Debye temperature (Qq is the volume per 
molecule), and V s i is the characteristic "modulation" of 
spin energy under long-wavelength acoustic deformation. 
At present, the physical mechanism of spin-lattice cou- 
pling is unclear, but the value of V s i — 2.6 K has been 
estimated from the available experimental data M . As a 
result, the estimate is (r s f r ) _1 ~ 2T • lO" 11 K (where T 
is the temperature in Kelvins). Here and below, we put 
h = ks = 1 and express all quantities, including relax- 
ation time, in the same units (Kelvins). Also, there is a 
contribution from Raman two-phonon processes, but at 
low temperatures the corresponding relaxation time t^ 
is very long: T$ r /r% ~ T 6 /(Mv 2 A§0 3 ) < 1 §§, and 
can be neglected. 

We also consider Orbach two-step relaxation via the 
excited levels S = 3/2 § (see Fig. |): 
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and for the parameters of V15 we have (t^ 7 ") ~ 
10 _8 exp (— Ai/T) K. Here, we assume that the spin- 
lattice matrix element V s i is of the same order as above 
(about 2.6 K). 

Along with triggering Orbach processes, the excita- 
tion of molecules to the level S = 3/2 leads to a vari- 
ation of the dipolar field exerted on a given molecule. 
As time goes by, some of the excited molecules relax 
back to S = 1/2, while other molecules go up to the 
level S = 3/2, and the dipolar field at a given point 
in the crystal fluctuates with time. In this paper, we 
use a mean-field approach to take into account the dipo- 
lar fields acting on molecules; it is justified since we are 
dealing with the case of relatively high temperatures (in 
comparison with the energy of dipolar interactions Tq) 
and long-range dipolar forces. Within the mean-field ap- 
proach, the dipolar field of the molecule with the spin 
S2 (equal to 3/2) can be imagined as a sum of the field 
created by a spin Si (equal to 1/2) and a field created by 
the spin S' = S2 — Si. Thus, the total dipolar field is a 
sum of two fields: the static demagnetizing field created 
by a uniform medium of spins 1/2, and a purely fluctuat- 
ing field h created by the spins S'. The spins S' are dis- 
tributed approximately uniformly over the sample at any 
instant, and their number iVi is small in comparison with 
the total number N of molecules, Ni — TVexp (— Ai/T), 
so the fluctuating field h at any instant obeys the Cauchy 
(Lorentz) distribution (Chapter IV, Ref. J8J): 



P(h) = l 
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with T = T (Ni/N) = r cxp(-Ai/T), where T 
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10~ 4 K is of order of the dipole-dipole interaction en- 
ergy in the ground state. Note that the fluctuating field 
h is measured against the total static field, including the 
static dipolar field. A comparison with Eqs. (Q) and (g) 
shows that at T > 0.2 K the distribution width T is much 
larger than 1/t s i, so that the fluctuating field h destroys 
coherence much faster than phonons do. Therefore, the 
fluctuating dipole-dipole field constitutes an important 
decoherence factor. 

To estimate the correlation time of the dipolar field 
fluctuations, we note that the field changes when excited 
molecules relax back to the S = 1/2 level (and, according 
to the principle of detailed balance, the same number of 
molecules go to the level S = 3/2). The transition from 
S = 3/2 to S — 1/2 proceeds via emission of phonons 
of energy Ai. The rate of this transition is proportional 
to A 3 (the number of phonons with energy Ai), and can 
be calculated using the Fermi's golden rule (Chapter 10, 
Ref. §): 

Mv 2 \0 J ' ( ' 

which satisfies the condition of detailed balance between 
the levels S = 3/2 to S = 1/2 (cf. Eqs. (|) and (§) repre- 
senting the rates of transitions "up" and "down" ) , so the 
level populations remain constant in time. During the 
time r c , a majority of the molecules situated initially in 
the state S = 3/2 relax to S = 1/2, and other molecules 
are excited, causing the field to fluctuate. Thus, r c is the 
correlation time for the fluctuating dipolar field. The es- 
timate gives r" 1 ~ 10~ 8 K for Vj.5, so that Tt c <C 1 at 
T < 0.5 K, i.e. the field fluctuations are fast in compari- 
son with their amplitude. 

Now, let us consider the hyperfine fields which consti- 
tute an important source of decoherence J7J. A typical 
time for fluctuations of the hyperfine field t„ is of order of 
the linewidth of nuclear magnetic resonance and can be 
estimated as dipole-dipole interactions between different 
nuclei J8|: l/r„ ~ (/i„//i e ) 2 ro, where n n , \x e are nuclear 
and electronic magnetic moments, respectively. There- 
fore, for the temperatures T > Ai/[2 In (/i„//i e )] ~ 0.2 K 
one has r n T 3> 1 and hyperfine fields can be considered as 
static for time intervals of order T _1 . As will be shown 
below, r defines the relaxation (decoherence) time, so 
that hyperfine fields can be combined with static demag- 
netizing fields to give some total static mean-field bias 
ho of the doublet levels. This bias is determined mainly 
by the hyperfine field exerted on a molecule, which is 
about Thf — 5 • 10~ 2 K [H (demagnetizing fields are 
weaker), and is of order of the tunneling splitting Ao. 
Therefore, for a large fraction of the molecules the lev- 
els \S Z = +1/2) and \S Z = —1/2) are rather close to 
resonance. This is radically different from the case of 
strongly anisotropic molecular magnets (such as Mni2 or 
Fes) where the ground-state tunneling splitting is much 
smaller than hyperfine fields. 



Finally, we consider the static dipolar interaction Tq ~ 
10~ 4 K between the molecules situated in the lowest four 
states (with S = 1/2). Longitudinal dipolar interac- 
tions (the terms SlS z ) are included in the mean field 
along with the static hyperfine field Thf , and can be ne- 
glected in comparison with the latter (since Tq <C Thf)- 
The terms S\S X etc. within the mean- field approxima- 
tion just change the tunneling splitting negligibly (since 
r <C Ao). But the flip-flop terms (S x Sy etc.) cannot 
be incorporated into the mean-field scheme. Flip-flop be- 
tween two molecules is a transition from the state \Sl = 
+1/2, S* = -1/2) to the state \S X Z = -1/2, S* = +1/2). 
The matrix element of this transition is of order Tq, but 
the energy difference between the initial and final states 
is determined by the difference in local mean fields act- 
ing on the two molecules, which is of order Thf 3> Tq. In 
this situation, known as Anderson localization, the levels 
of the molecule do not widen at all, and no spin diffu- 
sion is present. The localization can be lifted due to the 
dynamic change of the hyperfine field at the molecule, 
but this happens on a timescale t ~ r n . At tempera- 
tures T > 0.2 K the coherence is already lost at these 
times, due to thermoactivated dipolar field fluctuations 
(r-r n 3> 1). At lower temperatures, the mean-field ap- 
proach is not valid, and the intermolecular correlations 
should be taken into account. 

Summarizing the discussion above, the dipolar dy- 
namic fluctuations constitute an important source of de- 
coherence at 0.2 K< T < 0.5 K. Let us formulate now 
a model for magnetic relaxation under the fluctuating 
dipolar field h = h x e x + h y e y + h z e z . We consider a 
two-level system (the levels S z = ±1/2 for V15) with the 
static tunneling splitting Ao and static mean-field bias 
ho (the latter is governed mainly by the hyperfine static 
fields, since the demagnetizing fields are much weaker). 
The system is described by the density matrix p written 
in the basis formed by the levels S z = ±1/2. Its equation 
of motion is 



p = i[p, H] 



(5) 



where TC = — (Ao + h x )a x — h v a y — (ho + h z )a z is the 
Hamiltonian of the system (<J Xt y tZ are the Pauli's matri- 
ces). It can be conveniently written as 



x = -h y y - (A + h x )z 

y = h y x + (ho + h z )z 

z = (A + h x )x - (ho + h z )y 



(6) 



by introducing the variables: x = (pn — p 22 )/2, y = 
(p\2 + P2i)/2 and z — (p 12 — pii)/(2i). The static fields 
Ao and ho can be eliminated by two rotations of the co- 
ordinate frame: 

x = X cos tp — (Y cos Et + Z sin Et) sin tp (7) 

y = X sin ip + (Y cos Et + Z sin Et) cos tp 
z = -Y sin Et + Z cos Et, 



where s'mip = Aq/E, cosip = ho/E, E — ^/Aq + hg, and 
Eqs. (o) take the form 



V2X = (h 2a - h 3b ) Y - (h 2b + h 3a ) Z 
V2Y = sf2hxZ + (h 3b - h 2a ) X 
V2Z= (h 2b + h 3a ) X - V2hiY. 



(8) 



The random fields acting on the system are hi — 
h z cos (p + h x sin tp, h 2aM = y/2 ft.2,3 sin Et, and h 2by3b = 
v2 h 23 cos Et, where h 2 = — h z smtp + h x cost/? and h 3 — 
h y . As we discussed above, h x ,y,z ar e independent ran- 
dom fields, at any instant distributed with the same law 
(||); the same is true for /ii,2,3- Since Et c > Aqt c ^> 1, 
one can consider h 2a ,3a and h 2b ^ b as independent, and 
in Eq. (ph we have several independent fluctuating fields 
with the same Cauchy distribution and with very short 
autocorrelation time t c . 

Eqs. (ph can be imagined as describing the evolution 
of a system with the Hamiltonian H = Hi + H 2 + H 3 
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where R = (X,Y,Z), Hi = V2(h 2a - h 3b )S l: H 2 = 

— \ji(h, 2b + h 3a )S 2 , H 3 = hiS 3 , and the noncommuting 
matrices Si 2 3 are 



Si 



The formal solution of Eqs. (|9j) can be represented in 
a path-integral-like form, by splitting the time interval 
(0, t) into N ^> 1 equal pieces of length e = t/N: 

R(t) = exp [-ieH(t N -i)] . . . exp [-ieH (0)]R(0) (11) 

where t n — ne. Each of the matrices H is proportional 
to the fluctuating fields ft.1,2,3, so if we choose e <C T^ 1 
the Trotter decomposition formula pi can be applied to 
each factor: 



exp(-ieJ2H k ) =l[cxp (-ieH k ) + 0(e 2 



(12) 



where k — 1,2,3. The correlation time of all the fields 
^1,2,3 is Tc, so Hk(t) and Hk(t + e) in Eq. ( |l2| ) are decor- 
related if e ^ t c - Choosing r c <C e -C T _1 , each term 
in the products ( pd| ) and ( ^2j ) can be averaged indepen- 
dently over different realizations of the random processes 
represented by the fields hi 2 3 thus giving the answer: 



(X(t)) = X(0) exp (-2v / 2ri) 
(Y(t)) = Y(0) exp [-(V2 + l)Tt] 
(Z(t)) = Z(0) exp [-(V2 + l)Tt] 



(13) 



These results, together with Eq.dj), represent an ex- 
act solution of the problem. The situation considered 



here is similar to that found in spin resonance, and the 
results can be conveniently expressed in corresponding 
terms. The dynamics of the density matrix elements is 
represented as a sum of two terms: damped oscillations 
with the frequency E (with transverse relaxation rate 
T^ 1 = (\/2 + l)r) and pure damping (with longitudinal 
damping rate T-j - = 2^/2T). The decoherence times T\p 
both are of order r _1 . This holds in spite of the small- 
ness of t c , due to the peculiar properties of the Cauchy 
distribution: for Gaussian fluctuations with variance a 2 
we would have much smaller relaxation rate ct 2 t c (mo- 
tional narrowing pi). On the other hand, if r c were very 
large then the dipolar field would be almost static, and 
the decoherence time for a single molecule would be de- 
termined by hyperfine fields, as it is for Mn 12 or Fe 8 J7). 

Nevertheless, for V15 the decoherence rate is still small 
enough: T/A ~ 2 • 1(T 7 at T = 0.5 K, i.e. the system 
tunnels about 5,000,000 times before the tunneling oscil- 
lations are wiped out by decoherence. We emphasize that 
each tunneling in V15 is not a single-spin event: it takes 
place between the two states of the whole molecule. It 
is a tunneling of an antiferromagnetic system with small 
uncompensated spin, and all 15 spins are involved. 

Summarizing, we considered possible sources of deco- 
herence in V15 molecules between the states S z — ±1/2 
of ground-state doublets. We found that in the temper- 
ature region 0.2-0.5 K the main source of decoherence 
is the fluctuating dipolar field created by the molecules, 
which are thermally activated to the higher S — 3/2 level. 
Based on an exactly solvable model, a rather low deco- 
herence rate is found: about 5,000,000 tunneling events 
occur before the coherence is destroyed. Such a low de- 
coherence rate is unusual for magnetic systems of meso- 
scopic size at these temperatures. 

Authors would like to thank W. Wernsdorfer, I. 
Chiorescu and B. Barbara for helpful discussions. This 
work was partially carried out at the Ames Laboratory, 
which is operated for the U. S. Department of Energy 
by Iowa State University under Contract No. W- 7405-82 
and was supported by the Director for Energy Research, 
Office of Basic Energy Sciences of the U. S. Department 
of Energy. This work was partially supported by Russian 
Foundation for Basic Research, grant 98-02-16219. 



[3] 

[4] 
[5] 
[6] 

[7] 

[8] 
[9] 

[10] 



J. Tejada and R. Ziolo, Phys. Rev. Lett. 76, 3830 (1996); 

L. Thomas, F. Lionti, R. Ballou, D. Gatteschi, R. Sessoli 

and B. Barbara, Nature 383, 145 (1996); and references 

therein. 

D. Gatteschi, L. Pardi, A. L. Barra, A. Muller, and J. 

Doring, Nature 354, 465 (1991); D. Gatteschi, L. Pardi, 

A. L. Barra, and A. Miiller, Molecular Engeneering 3, 157 

(1993). 

I. Chiorescu, W. Wernsdorfer, A. Muller, H. Bogge, and B. 

Barbara, submitted to Phys. Rev. Lett. 

P. W. Shor, SIAM J. Comput. 26, 1484 (1997); B. E. Kane, 

Nature 393, 133 (1998). 

D. Garanin and E. Chudnovsky, Phys. Rev. B 56, 11 102 

(1997); A. Fort, A. Rettori, J. Villain, D. Gatteschi, and 

R. Sessoli, Phys. Rev. Lett. 80, 612 (1998); F. Luis, J. 

Bartolome and J. Fernandez, Ph ys. Rev. B 57, 505 ( 1998); 



M. N. Leuenberger and D. Loss, cond-mat /9810156 



N. V. Prokof'ev and P. C. E. Stamp, |sond-mat/9511015| ; 
J. Low Temp. Phys. 104, 143 (1996); Phys. Rev. Lett. 80, 
5794 (1998). 

A. Abragam, The Principles of Nuclear Magnetism 
(Clarendon Press, Oxford, 1961). 

A. Abragam and B. Bleaney, Electron Paramagnetic Reso- 
nanse of Transition Ions (Clarendon Press, Oxford, 1970). 
L. S. Schulman, Techniques and Applications of Path Inte- 
gration (Wiley, New York, 1981) 



FIG. 1. Sketch of the V15 molecule. The spins B and B' of 
the upper hexagon form a very stiff antiferromagnetic dimer: 
Jbb' ~ 800 K. The spin A from the central triangle is coupled 
antiferromagnetically to both B and B' , so Jab and Jab' 
cancel each other to a large extent (Jab, Jab' "C Jbb')- This 
frustration effectively decouples the triangle spins from the 
hexagons. Due to C3,, symmetry of the molecule, the situation 
is the same for all the triangle spins. 
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